Quarter-filled supersolid and solid phases in the extended Bose-Hubbard model 
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We numerically study the ground state phase diagram of the two-dimensional hard-core Bose- 
Hubbard model with nearest (Vi) and next nearest neighbor (V2) repulsions. In particular, we focus 
on the quarter-filled phases where one supersolid and two solid phases are observed. Using both 
canonical and grand canonical quantum Monte Carlo (QMC) methods and a mean field calculation, 
we confirm the existence of a commensurate supersolid at quarter boson filling. The nature of 
the commensurate supersolid will be discussed. Only one kind of the supersolid phase is found 
energetically stable despite of the two possible diagonal long range orderings for the solid phase. 
The competition between the two solid phases manifests as a first order phase transition around 
2V2 ~ Vi. The change of order parameters as functions of the chemical function is also presented. 

PACS numbers: 75.10.Jm, 75.45.+J, 05.30.Jp, 75.40.Mg 



Bose-Hubbard model has recently attracted a lot of 
attentions for the possibility of observing the supersolid 
phase Jl[ either in optical lattices or in magnetic sys- 
tems !3|. The simultaneous breaking of both transla- 
tional and 17(1) gauge symmetry is a delicate state of 
matter that only recently experimental evidence of pos- 
sible supersolid phase is provided by the measurement 
of the non-classical rotational inertia (NCRI) 4] in ro- 
tating solid 4 He However, the observed NCRI may 
be attributed to the superfLow in between microcrystal 
interfaces and the issue is still largely controversial Q. 
Thanks to the technological advance in trapping atoms 
or even molecules at very low temperature in optical lat- 
tices, it provides an ideal testing ground for the searching 
of the supersolid phase. Hubbard model of hard-core bo- 
son in the frustrated triangular lattice Q and soft-core 
boson 0] in square lattice are among the possible can- 
didates. Besides the optical lattice experiments, various 
magnetic systems has been suggested [3] to be candidates 
of the realization of spin supersolid in some carefully cho- 
sen parameter regime. In the case of spin S = h, where 
the system is equivalent to the hard-core Bose-Hubbard 
model, the spin supersolid represents the state of spatial 
modulation of the in-plane spin projection. 

Hard-core Bose-Hubbard model with only nearest 
neighbor (nn) interaction on a square lattice has ground 
states of supcrfluid ordering and checkboard solid order- 
ing. On the other hand, sufficiently large next nearest 
neighbor (nnn) interaction leads the checkerboard or- 
dering replaced by a striped ordering which can coex- 
ist with a superfiuid to form a striped supersolid around 
half- filling 0. Recently we revisited the model [l(| and 
found a new solid phase with 'star' ordering (Fig. QJ,) 
at quarter-filling. Interestingly, this star solid (A-phase) 
ordering can also coexist with the superfiuid ordering 
to form a star supersolid phase at and around quarter- 
filling. It has also been suggested ll| that another solid 
(B-phase) order can be stabilized when 2V2 > Vi. How- 
ever, whether the supersolid ground state with this solid 



ordering is stable is unclear and the complete ground 
state phase diagram including the B-phase solid is still 
absent. These are the questions we attempt to address 
in this work. Furthermore, there are concerns about the 
star supersolid phase at exactly quarter-filling for it con- 
tradicts the belief that the existence of commensurate 
supersolid 
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Il2j is impossible. We will show in this 
paper that the quarter-filled star supersolid is commen- 
surate but agrees with the notion of vacancy supersolid. 

Specifically, we study the extended Bose-Hubbard 
model on a 2D square lattice with the Hamiltonian 



H b = -tyi(blbj+bjb])+Vi } j n i nj-\-V2 } j njnj-fiyn 

i,j i,j i,j i 

(1) 

where bQr) is the boson destruction (creation) operator 
and j2 nn (J2 nnn ) sums over the (next) nearest neighbor- 
ing sites. For convenience we fix the energy scale by set- 
ting t = 1 throughout this paper. By taking the transfor- 
mation — > and n —> S z + \, this Hamiltonian can be 
mapped to a spin XXZ model with nn and nnn exchange 
couplings under a magnetic field h = /1 — | (Vi + V2) (z 
is the coordination number). At half- filling, the ground 
state of the Hamiltonian Hb is a checkerboard solid (char- 
acterized by wave vector Q=(7r, it)) for strong nn cou- 
pling V\ , or a striped solid (characterized by wave vector 
Q=(7r,0) or (0,7r)) for strong nnn coupling V2 [§j]. For 
V\ ~ 2V2, quantum frustration destabilizes both solid 
orders and leads to a uniform superfiuid phase. No su- 
persolid phase is found stable at half-filling. Away from 
half-filling, however, a striped supersolid is observed for 
dominating nnn interaction 2V2 > Vi, while the checker- 
board supersolid phase is still unstable against phase sep- 
aration for large nn interaction Vi > 2V2. It has been ar- 
gued that the motion of domain walls reduces the ground 
state energy so that the checkerboard supersolid is ener- 
getically unstable [§]. 

Further away from half-filling, new types of solid in 
the vicinity of quarter-filling, as well as three-quarter- 
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filling, are found very recently 13, HI, E3] as mentioned 
above. The A-phase quarter-filled solid (see Fig. [1]) has 



finite structure factor S(Q)/N = J2i 



)/N 2 at 



wavevectors Q = (tt, tt), (tt, 0) and (0, tt), which implies 
a star-like occupation pattern. Moreover, doping the star 
solid with extra bosons yields a star supersolid via a sec- 
ond order phase transition. The formation of domain 
walls is no longer energetically favorable and hence the 
star supersolid is stable upon doping instead of phase sep- 
aration. More interestingly, this star supersolid persists 
even at exact quarter-filling for wide range of parameters 
V\ and Vi (see Fig . [8]). This result seems to contradict 
a recent proof [l5j that the necessary condition for su- 
persolidity is to have zero-point vacancies or defects and 
no commensurate supersolid is possible. Ref. 11, [l2j] 



claim they do not observe a quarter-filled star supersolid 
in their QMC data. Unfortunately, no data of the struc- 
ture factor and superfluidity are presented in the right 
Vi , V2 parameter regime where we observed the quarter- 
filled supersolid. One natural question arises about the 
discrepancy is whether a canonical approach, used in ref. 
[13] where particle number is fixed, leads to different re- 
sults. In order to clarify the issue, we provide further 
evidence using Green's Function Monte Carlo (GFMC) 
to support the existence of quarter-filled supersolid. 



STOCHASTIC SERIES EXPANSION 

In order to clarify the issue of the presence of the 
SS phase at n — 0.25, we numerically study the model 
with both the grand canonical and canonical approaches. 
The grand canonical calculation is carried out using the 
standard stochastic series expansion (SSE) Monte Carlo 
method implemented with directed loop algorithm [l3j |. 
While SSE works on grand canonical ensemble, in order 
to fix the density n we scan the chemical potential fi to 
find an average density (n) = 0.25 with an uncertainty 
less than 0.01. It turns out, as shown below, this method 
generates the same result as the GFMC. In SSE, the su- 
perfluidity, given by p x{y) = (W 2 , y) )/A(3t, is computed by 
measuring the winding number fluctuation as usual. 



GREEN'S FUNCTION MONTE CARLO 

The GFMC starts with a Jastrow variational wave 
function which is defined by applying the density Jas- 
trow factor to a state with all the bosons condensed into 
the q — state 



= exp 



(2) 
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(a) Solid A 



(b) Solid B 



where |$ ) = (EA^I ) is the non-interacting boson 
ground state with iV particles, and v^j are parameters 



FIG. 1: Classical lattice structures of star solid (a) A and (b) 
B phase. 



that can be optimized to minimize the variational en- 
ergy In order to take the hard-core constraint into 
account, configurations with more than one boson on 
a single site are projected out from Eq. [2l The wave 
function \^>j) in Eq. [2] was shown to be able to turn a 
non-interacting bosonic state into a Mott insulator if a 
long-range Jastrow factor is included. In our recent work, 
we have also shown that the supersolid and solid phases 
can also be described in the same wave function. How- 
ever, the number of the variational parameters in Vij 
grows exponentially with the lattice size which costs a 
lot of computation time for an optimized wave function. 
Instead of including parameters of all range, we use a 
Gutzwiller projection factor exp {X^eAff 77- *} ^° enhance 
the diagonal order. Here g is a variational parameter 
which controls the diagonal order specified by the sub- 
lattice A. Obviously, this factor can stabilize the solid 
and gives reasonably good trial energy when V\ and V2 
are large. We found that low variational energy can be 
acquired for the A-type supersoild and solid without the 
projection factor. On the other hand, relative stable wave 
function is found by including the projection factor for 
the B-type solid. 

In order to investigate the exact ground state proper- 
ties, Green's function Monte Carlo method is employed to 
improve the variational results. In this work, we choose 
multi-walkers stochastic reconfiguration method to pre- 
vent the simulation from blow-up or dead-ends in the 
large power limit. To benchmark our method we compare 
our GFMC data (points) with the exact results (lines) 
without using Monte Carlo for an 6 x 6 lattice in Fig. [5] 
We can see that the GFMC results are consistent with the 
exact ones. The underestimation of the diagonal order in 
the variational wave function is corrected as the iteration 
increases. We have verified that the same ground state 
properties can also be obtained with the wave function 
with diagonal order (g 7^ 0). 

Now we present the GFMC results for larger lattice. 
In Fig. [3] we show the structure factors and condensate 
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FIG. 2: GFMC results (solid symbols) of (a) energy, the 
structure factor of wavevectors (b) (tt, 0) and (c) (n,ir), and 
(d) condensate density iVfc=o (square) as functions of iteration. 
Solid lines are exact results without using Monte Carlo. Here 
Vi = 2, V2 = 4.5, lattice size is 6x6 and n = 0.25. The guiding 
function are the same as trial wave function with 1)1,0 = 0.9, 
»i 1 = 0.83. 4000 walkers are used in the GFMC calculation. 



Nk=o = b\. =0 bk=o as a function of iteration. The trial 
wave function with A-type diagonal order is used as the 
trial wave function. In order to check that the ground 
state can be reached regardless of the choice of g, we 
present results of g = 0.2 and g = 0.4. The wave function 
is optimized for all parameters before GFMC is applied. 
The corresponding data obtained by SSE (dash lines) are 
also shown for comparison. As we can see that although 
the optimized wave function overestimated the structure 
factor and underestimated the condensate density, con- 
sistent results with SSE are obtained as the iteration in- 
creases. Fig. 2]are similar calculations for V2=3.4 and 4.4 
corresponding to superfluid and A-type solid phases, re- 
spectively. The data clearly show the convergence of the 
GFMC approach and consistent with the SSE results. 
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FIG. 3: (Color online) (a) The GFMC results of the structure 
factor of wavevectors (n, 0), (71", n), and (b) condensate density 
Nk=o (square) as functions of iteration. Dash lines are SSE 
results. Vi = 4, V2 = 3.8, lattice size is 12 x 12 and n = 0.25. 
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FIG. 4: (Color online) (a) The GFMC results of the structure 
factor of wavevectors (n, 0), (n, n), and (b) condensate density 
Nk=o (square) as functions of iteration. Dash lines are SSE 
results. Here Vi = 4, V2=3A and 4.4. Lattice size is 12 x 12 
and n — 0.25. 



MEAN FIELD THEORY 

To further investigate the effect of quantum fluctua- 
tion of the model, we also obtain the ground state phase 
diagram using a simple mean field approach. A mean 
field wavefunction 

|*)mf = nW 1 >v*+^l >v*) ( 3 ) 

is given to represent the superfluid, star solid and star 
supersolid phases. Here a denotes the index of a 2x2 
unit cell while i = 1, 2, 3, 4 is the sublattice index inside 



the cell. \l)i t a (|0)j, a ) is the occupied (empty) Fock state 
at i-th site of the unit cell a while ut{vi) is the corre- 
sponding variational parameter. The energy of \^)mf is 
then minimized to obtain the ground state wavefunction. 
This MF wavefunction successfully predicts the existence 
of quarter-filled star supersolid as well as the superfluid 
and star solid as shown in the inset of FigO 



QUARTER-FILLED SUPERSOLID 

FigG3 shows the result of SSE and GFMC at Vi = 4 
and V\ = 8 for small lattice size. The agreement be- 
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FIG. 5: (color online) Comparison of data obtained from the 
SSE (filled symbols) and GFMC (open symbols) for lattice 
size L = 12. Note that GFMC measures condensate density 
Nk=o while SSE measures superfluidity p for the signature of 
superfluid phase. Temperature used in the SSE is T/t=0.02. 
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FIG. 6: (color online) Finite size analysis of the order param- 
eters p x and S(iv, 0) from SSE for Vi = 2.0. 



tween both approaches is clear and verifies that our at- 
tempt to fix the boson density n to 0.25 in the grand 
canonical SSE does not lead to any measurable discrep- 
ancy of the physical quantities we are interested in with 
the canonical GFMC. The coexistence of superfluid or- 
der and star crystal structure for 3.5 < V2 < 4.0 in FigO 
(a) clearly signals the supersolid phase at quarter-filling. 
The uniform superfluid develops spatial modulation, i.e. 
becomes a supersolid, continuously as V2 increases and 
gradually loses its superfluidity at the same time until 
it eventually becomes a star solid. To demonstrate that 
the supersolid ground state survives at larger lattice size 
and at thermodynamic limit, a finite size analysis is car- 
ried out in Fig. [6] for Vi = 2. S(ir, 0) scales to zero 
in the superfluid phase whereas p x scales to zero in the 
solid phase. Only in the supersolid phase (V2 — 4.4 in 
Fig. [6]) that both S(tt, 0) and p x scales to finite values. 
It is remarkable that simple MF theory also correctly 
predicts the existence of quarter-filled SS and quantum 
fluctuation does not destroy the long-range orders in the 
SS phase, in contrast to the case of Kagome lattice [lj|. 
Except that the solid A and B phase are indistinguish- 
able in the MF level, MF theory successfully reproduce 
all phases at quarter-filling as shown in the inset of Fig. 

m 

In the first glance, the existence of a supersolid at com- 
mensurate density contradicts the notion of vacancy su- 
persolid. Prokof 'ev and Svistunov 15j have proved that 
superfluidity has zero probability to occur in commen- 



surate solids in nature, or in other words, the necessary 
condition for the supersolidity is the present of vacancies 
or defects. This is due to the asymmetry between vacan- 
cies and interstitials. However, this result, as admitted 
in their paper, does not apply to systems with explicitly 
broken translation symmetry, such as in lattice models 
where commensur ability is enhanced by hand. In our 
system, the commensurability of quarter-filled supersolid 
can always be ensured by adjusting the chemical poten- 
tial provided there is no phase separation. Vacancies and 
interstitials (means bosons at the sublattice sites) arise 
from quantum fluctuations do not form bound pairs and 
are free to move away from each other that, therefore, 
leads to superflow. The interstitial-vacancy symmetry is 
generally absent in nature but is preserved in this case 
by the external potential that fix bosons to locate only 
at the lattice points. Based on the measured structure 
factors S(tt, 7t), S(tt, 0) , S(0, it) and boson density n, 
one can easily deduce the boson densities on each sublat- 
tices. Fig|7] shows that bosons do not localized at only 
one sublattice but have finite occupation in all four sub- 
lattices. The data indicates that there are more than 30 
percent of vacancies in the quarter-filled supersolid and 
in consistent with the picture of vacancy supersolid. 

To gain more insight of the quarter-filled supersolid 
phase, we plot in Fig[7]the order parameters as functions 
of the chemical potential fj,. Starting from small density, 
the ground state is an uniform superfluid in which all 
sublattice densities are equal. While increasing [i until 
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FIG. 7: (color online) SSE result of the (a) average density 
n, (b) sublattice densities, (c) superfluidity and (d) structure 
factors as functions of chemical potential for Vi = 2.0 and 
V2 = 4.6. Open (filled) symbols denote the result for a lattice 
of size 24x24 (48x48). 



quarter-filled, the system undergoes a second order phase 
transition to a star SS. Note that although there is a small 
dip of superfluidity in the SS state, it does not reduce to 
zero even at n = 0.25. There is no any indication that 
n = 0.25 is a special density that acquires particular 
treatment like the canonical calculation. Furthermore, 
no noticeable change is observed when doubling the lat- 
tice size to 48x48 (filled symbols in FigJT]). Doping more 
bosons will destabilize the star SS phase because of the 
nnn repulsion that leads to a striped SS via a discontin- 
uous phase transition between the two different broken 
translation symmetries. Striped SS is also observable at 
n = 0.25 when V\ is small enough. On the other hand, 
while there are two kinds of star solid (A and B phase) , 
a natural question is whether SS of both kinds exist at 
or away from quarter-filled . We will address this issue 
in the following section. 
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FIG. 8: Phase diagram of n = 0.25. The lines are guides to 
the eyes. Solid (dotted) lines represent second (first) order 
phase transition boundaries. Black circles (squares) denote 
the data obtained from the SSE (GFMC) for L=28 (L=16). 
The inset shows the MF result. It is noted that star solid A 
and B phase cannot be distinguished in the MF level. 

SUPERSOLID A AND B PHASE 

It is interesting to note that at quarter-filling, the clas- 
sical ground state of the frustrated Hamiltonian H is 
highly degenerate. Translating any lines of bosons in Fig. 
QJ, by one lattice constant of the solid A phase along, say, 
x direction will create a domain wall with no energy cost. 
When translating alternate lines of bosons generate a B 
phase solid. Enormous number of ways to create domain 
walls implies that the classical ground state has macro- 
scopic degeneracy. This degeneracy, however, is lifted 
by quantum fluctuation that yields a ground state of A 
phase if V\ < 2V2 or B phase otherwise. This is another 
typical example of the order by disorder phenomenon Q. 
As discussed in previous sections, it leads to a star SS A 
phase in a wide parameter range, similar to the scenario 
of the one-third filling SS in the frustrated triangular lat- 
tice [7j. On the other hand, the existence of star SS B 
phase has not been clarified. Here we complete the phase 
diagram for larger V\ where solid B phase can be stabi- 
lized at quarter-filling. Fig. [5] shows that, for V\ > 6, 
when increasing from small nnn interaction V2, the SF 
phase changes to the solid B phase via a first order phase 
transition without passing an intermediate SS B phase. 
Star solid A phase emerges when further increasing V2 
until 2V2 > Vi . Furthermore we do not observe any SS B 
phase away from n = 0.25 as shown in Fig. [9] at a repre- 
sentative Vi = 8.0 and V2 — 3.5. Instead there is a first 
order phase transition from the gapped solid B phase to 
a uniform SF phase and implies that the SS of B phase 
symmetry is unstable towards phase separation. This can 
be understood by the simple argument of domain forma- 
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FIG. 9: (Color online) SSE result for Vi = 8.0 and V 2 = 3.5. 
The ground state at n = 0.25 is a B phase solid. Doping 
away from n = 0.25 leads to a first order phase transition to 
SF phase. No B phase supersolid is found. 

tion as discussed in the case of hardcore checkerboard 
SS 0] . Addition holes (bosons) in the quarter-filled solid 
phase B tend to line up to form a domain wall in which 
the gain in kinetic energy is larger than that of isolated 
holes (bosons). The SS B phase is therefore unstable 
towards phase separation. 

SUMMARY 

We have presented a comprehensive numerical study 
on the extend hardcore Hubbard model, in particular, 
on the SS phases at the quarter-filled density. Based on 
results of the SSE, GFMC and MF calculations, we pro- 
vide convincing evidence for the existence of star SS A 
phase at exact n = 0.25, in contrasts to previous study. 
We clearly show that the SS A phase found at n = 0.25 is 
consistent to the notion of vacancy SS. The star SS phase 
is a consequence of the order by disorder phenomenon by 
which the ground state degeneracy is lifted. Further- 
more, the physical natures the solid A and B phase are 
also studied. Although classically degenerate, these two 
states compete with each other through quantum fluctu- 



ations and the final stability depends on the competition 
between nn (Vi) and nnn (V2) interactions. We present a 
complete phase diagram including also the solid B phase. 
The SS B phase on the other hand is found unstable to- 
wards phase separation due to the kinetic energy gain of 
domain formation. 
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